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Abstract 

Synchronization of coupled continuous-time linear systems is studied 
in a general setting. For identical neutrally-stable linear systems that 
are detectable from their outputs, it is shown that a linear output feed- 
back law exists under which the coupled systems globally asymptotically 
synchronize under all fixed (directed) connected network topologies. An 
algorithm is provided to compute one such feedback law based on indi- 
vidual system parameters. The dual case, where individual systems are 
neutrally stable and stabilizable from their inputs, is also considered and 
parallel results are established. 

1 Introduction 

In [20] we have shown, for identical discrete-time linear systems that are de- 
tectable (stabilizable) from their outputs (inputs) and neutrally stable, that a 
linear feedback law exists under which the coupled systems globally asymptot- 
ically synchronize for all fixed (directed) connected network topologies. There 
we have also provided an algorithm to compute such feedback law based on 
individual system parameters. In this companion paper we provide counterpart 
results for continuous-time linear systems. 

1.1 Background 

"The main issue in studying the synchronization of coupled dynamical systems 
is the stability of synchronization. As in all cases where stability is the issue, the 
question whose answer is sought is Under what conditions will the individual 
systems synchronize? In a simplified yet widely-studied scenario, where the 
individual system dynamics are identical and the coupling between them is 
linear, studies focus on two ingredients: the dynamics of an individual system 
and the network topology. Starting with the agreement algorithm in 19J a 
number of contributions [8j [10l HH [TJ [12] have gathered around the case where 
the weakest possible assumptions are made on the network topology at the 



expense of restrictive individual system dynamics. It was established in those 
works on multi-agent systems that when the individual system is taken to 
be an integrator and the coupling is of full-state, synchronization (consensus) 
results for time-varying interconnections whose union^] over an interval are 
assumed to be connected instead of that each interconnection at every instant 
is connected. 

"Another school of research investigates networks with more complicated 
(nonlinear) individual system dynamics [lHld^]- When that is the case, the 
restrictions on the network topology have to be made stricter in order to ensure 
stability of synchronization. Generally speaking, more than mere connectedness 
of the network has been needed: coupling strength is required to be larger than 
some threshold and sometimes a symmetry or balancedness assumption is made 
on the connection graph. Different (though related) approaches have provided 
different insights over the years. The primary of such approaches is based on the 
calculations of the eigenvalues of the connection matrix and a parameter (e.g. 
the maximal Lyapunov exponent) depending on the individual system dynamics 
pMl IT51 [B]. In endeavor to better understand synchronization stability tools 
from systems theory such as Lyapunov functions [HE], passivity [Ml [3l H7l [23] . 
contraction theory [16 , and incremental input-to-state stability (<5ISS) theory 
[5] have also proved useful. "0 

1.2 Contribution 

In this paper we study two dual problems. In the first case we consider the 
following individual system 

Xi = Axi , yi = Cxi , (1) 

where A is assumed to be neutrally stable and pair (C, ^4) detectable, and design 
a linear output feedback gain L that synchronizes any fixed connected network 
of any number of coupled replicas of {!]) . Such L guarantees the synchronization 
of p individual systems when coupled as 

p 

Xi = Ax t + Hi (Vj -Vi)- 

3 = 1 

As the dual problem we consider 

±i = Axi + , (2) 

where A is assumed to be neutrally stable and pair (^4, B) stabilizable, and design 
a linear feedback gain K that synchronizes any fixed connected network of any 

1 By union of interconnections we actually mean the union of the graphs representing the 
interconnections. 

2 Borrowed from [20] . 
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number of coupled replicas of @ . Such K guarantees the synchronization of p 
individual systems when coupled as 



To the best of our knowledge, feedback design (in such a general setting) in order 
to guarantee synchronization under arbitrary (fixed) interconnections is a nov- 
elty of our work. It is worth noting that our main theorems make a compromise 
result between the two previously mentioned cases (i) where synchronization is 
established for very primitive individual system dynamics, such as that of an 
integrator, but under the weakest conditions on the network topology and (ii) 
where the network topology has to satisfy stronger conditions, such as that the 
coupling strength should be above a threshold, for want of achieving synchro- 
nization for nonlinear individual system dynamics. 

1.3 Organization 

The remainder of the paper is organized as follows. In the next section we 
provide notation and some preliminaries. Then we formally state our problems 
in Sections [3] and |U Section O is where we establish our key result which we will 
later use to solve the problems we aim at. In Section [5] we provide an algorithm 
to design output feedback gain that we seek for synchronization and prove that 
it works. Then, in Section [71 we design a state feedback gain that solves the 
dual problem. 

2 Notation and definitions 

Let N denote the set of nonnegative integers and K>o set of nonnegative real 
numbers. Let | • | denote 2-norm. Identity matrix in R nx ™ is denoted by I n . A 
matrix A G M. n x ™ is Hurwitz if all of its eigenvalues have strictly negative real 
parts. A matrix S G M" x ™ is skew- symmetric if S + S T = 0. Given C G R mx " 



and A G R nXn , pair (C, A) is observable if [C T A T C T A 2T C T . . . A^ n ^ T C T ] 



is full row rank. Pair (C, A) is detectable (in the continuous-time sense) if 
that Ce At x = for some x G W 1 and for all t > implies lim^oo e At x = 0. 
Given B G R nxm and A G W IX '\ pair (A, B) is controllable (stabilizable) if 
(B T , A T ) is observable (detectable). Matrix A G W ixn is neutrally stable (in 
the continuous-time sense) if it has no eigenvalue with positive real part and 
the Jordan block corresponding to any eigenvalue on the imaginary axis is of 
size oneH Let 1 G R p denote the vector with all entries equal to one. 

3 Note that A is neutrally stable iff there exists a symmetric positive definite matrix P such 
that A T P + PA < 0, 0. 
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Kronecker product of A £ R mx ™ and B £ W K( ' is 



ax\B 



a m \B 



a ln B 



iB 



Kronecker product comes with the properties (A (g> B)(C ® D) — (AC) <E) (BD) 
(provided that products AC and BD are allowed) A®B + A®C = A®(B + C) 
(for B and C that are of same size) and (A <g> B) T = A T (8 ~B T . 

A (directed) graph is a pair (A/", *4) where A/" is a nonempty finite set (of 
nodes) and A is a finite collection of pairs (arcs) (rii, nj) with rii, rij £ AT. A 
path from ni to ri£ is a sequence of nodes {n\, n2, ■ ■ ■ , ^} such that (rii, ^i+i) 
is an arc fori G {1, 2, £—1}. A graph is connected if it has a node to which 
there exists a path from every other node0 

The graph of a matrix T := [7^] G W pxp is the pair (A/", .4) where J\f = 
{n\, n2, . . . , Up] and (n^, n^) £ A iff 7jj > 0. Matrix T is said to be connected 
(in the continuous-time sense) if it satisfies: 

(i) 7y > for i ^ j; 

(ii) each row sum equals 0; 

(iii) its graph is connected. 

For connected T, it follows from definition that A — is an eigenvalue with 
eigenvector 1 (i.e. Tl = 0.) Moreover, all the other eigenvalues have real 
parts strictly negative. Let r T be the left eigenvector of eigenvalue A = (i.e. 
r T T = 0) with r T l = 1 . Then lim^oo e r * = lr T . 

Given maps & : IR>o — ► M. n for i = 1, 2, . . . , p and a map £ : M>o — > K™, 
the elements of the set {&(■) : i = 1, 2, . . . , p} are said to synchronize to £(•) if 
~ ~> as f -> 00 for all i. 



3 Problem I 

We now formalize our first problem. 



3.1 Systems under study 

We consider p identical linear systems 

±i = Axi + Ui , y i = Cx ll i = l,2,...,p (3) 

where Xi £ K n is the state, U{ £ R" is the input, and yt £ M. m is the output of 
the ith system. Matrices A and C are of proper dimensions. The solution of 

4 Note that this definition of connectedness for directed graphs is weaker than strong con- 
nectivity and stronger than weak connectivity. 
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ith system at time t > is denoted by x< (t) . In this paper we consider the case 
where at each time instant only the following information 



p 



(4) 



is available to ith system to determine an input value where 7^ are the entries 
of the matrix F € K pxp describing the network topology. Nondiagonal entries of 
r are nonnegative and each row sums up to zero. That is, the coupling between 
systems is diffusive. 

3.2 Assumptions made 

We make the following assumptions on systems ([3]) which will henceforth hold. 

(Al) A is neutrally stable. 
(A2) (C, A) is detectable. 

3.3 Objectives 

Our first objective is to show that there exists a linear feedback law L € ]^« xm 
such that, for all p and connectedT £ R pxp , solutions of systems ^ with Ui = 
Lzi, where Zi is as in globally (i.e. for all initial conditions) synchronize to 
a bounded trajectory. Our second objective is to devise an algorithm to compute 
one such L. 

4 Problem II 

In this section we state the second problem, which, as noted earlier, is the dual 
of the first. 

4.1 Systems under study 

Consider p identical linear systems 



where Xj £ W 1 and 6 M. m . Matrices A and B are of proper dimensions. We 
consider the case where at each time instant the following information 



Xi = Axi + Bu. 



•l ; 



i = 1, 2, p 



(5) 



v 




(6) 



is available to ith. system to determine an input value. 
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4.2 Assumptions made 

We make the following assumptions on systems ([5]) which will henceforth hold. 

(Bl) A is neutrally stable. 
(B2) (A, B) is stabilizable. 

4.3 Objectives 

Our first objective regarding the dual problem is to show that there exists a 
linear feedback law K £ R mx ™ such that, for all p and connected T € W xp , 
solutions of systems ([5]) with Ui — Kz i7 where is as in globally (i.e. for 
all initial conditions) synchronize to a bounded trajectory. Our second objective 
is to devise an algorithm to compute one such K . 

5 A special case 

Before we attempt to solve Problems I and II, we first establish a preliminary 
result to be resorted later. Consider the following coupled systems 

p 

t .SC • Il'liy^-.J^ i = l,2,...,p (7) 

3=1 

where & £ R" is the state of the ith system, S e M. nxn , &nd H e M. mxn . We 
make the following assumptions on systems ([7]) which will henceforth hold. 

(CI) S is skew-symmetric. 
(C2) (H, S) is observable. 
(C3) T := [jij] is connected. 

Below we provide our first result. 

Theorem 1 Consider systems 0. Let r 6 M. p be such that r T T = and 
r T l = 1. Then solutions for i = 1, 2, . . . , p, synchronize to 

6(0)" 
6(0). 

Proof. Consider matrix T - lr T . Observe that (r - lr T ) k = T k + (-l) fe lr T 
for k € N. For <6l therefore we can write 

e^ 1 ^ = L p + t{T-lr T ) + ^{T-lr T Y + ... 

= (/ P + tr + ^ + ...)-(^-^ + ...) 

= e rt - (1 - e-*)lr T . 



:= (r T ® e st ) 
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Consequently lim^^ e^ r lrT ) 1 = 0. We deduce therefore that T — lr T is Hur- 
witz. Since T — lr T is Hurwitz, there exist symmetric positive definite matrices 
P,QeR pxp such that 

-Q=(T-lr T ) T P + P(T-lr T ). (8) 

Define positive semidefinite matrices P :— (I p — lr T ) T P(I p — lr T ) and Q := 
(Ip — lr T ) T Q(I p — lr T ). Now pre- and post-multiply equation (JSJ) by (I p — lr T ) T 
and (I p — lr T ), respectively. We obtain 

-Q = (I p -lr T f(T-lr T ) T P(I p -lr T ) 

+ {Ip-lr T ) T P{T-lr T )(Ip-lr T ) 
= T T P(I p - lr T ) + (Ip - lr T ) T PT 

= T T (I P - lr T ) T P(I p - lr T ) + (I p - lr T ) T P(I p - lr T )T 
= T T P + PT . 

We now stack the individual system states to obtain x := [£f £j • • ■ £j] T . We 
can then cast ([7]) into 

x = (I p ® S + T ® P T P)x . (9) 

Define V : R pn -> R> as F(x) := x T (P ® J„)x. Differentiating F(x(*)) with 
respect to time we obtain 

V(x) = yL T (I p ®S T + T T (g>H T H)(P(g>I n )x 

+x T (P ® i n )(i p ® 5 + r $ p t p)x 
= x T (P® (5 T + 5) + (r T p + Pr)® J ff T J ff)x 

= -x T (Q®P T P)x. (10) 

Thence V(x) < for both Q and H T H (and consequently their Kronecker 
product) are positive semidefinite. 

Given some ( S R p ™, let A" C W n be the closure of the set of all points r] 
such that r\ — (lr T (g) e St )( for some t > 0. Set A" is compact for it is closed by 
definition and bounded due to that £ is fixed and S is a neutrally-stable matrix. 
Having defined X, we now define 

fl := { V e R pn : (lr T ® I n )r) £ X , V(rj) < V(()} ■ 

Let us show that f2 is forward invariant. Observe that 

j t ((lr T ® J„)x(t)) - (lr T ®I n )(Ip®S + T® H T H)x(t) 

= (lr T ® S + lr T T ® H T H)x(t) 

= (lr T ® S')x(t) 

= (/p®S')(lr T ®/„)x(t). 
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We therefore have 



(lr T ® Z»)x(t) = (lr T ® e st )x(0) (11) 

which in turn implies that if (lr T (gi 7„)x(0) G A" then (lr T ® Z„)x(i) £ A" for 
all t > 0. Likewise, if V(x(0)) < F(C) then V(x(i)) < F(C) for all * > thanks 
to flO]). As a result, if x(0) G Jl then x(i) G f2 for all t > 0, that is, £1 is forward 
invariant with respect to (J9j> . 

Set il is closed by construction. To show that it is compact therefore all we 
need to do is to establish its boundedness. Let 

a := sup 1 77 - (lr T <g Z„)7?| . 

V(77)<V(C) 

If we go back to the definition of V we immediately see that a < 00. Now let 

b :— sup \u>\ . 

Since X is bounded, b < 00 as well. Now, given any 77 G SI we have 1 77 — (lr T <g) 
Zn)?7| < a- Hence we can write 

|r?l < a+ |(lr T ®/„)?7| 
< a + sup |w| 

= a + b. 



Therefore il is bounded. Having shown that il is forward invariant and compact, 
we can now invoke LaSalle's invariance principle [9l Thm. 3.4] and claim that 
any solution starting in il approaches to the largest invariant set W C {?/ G £1 : 

v(v) = o}. 

Let now rj(-) be a solution of © such that rj(t) G W for all t > 0. Given 
some r > 0, since V r (?y(T)) = 0, we can write 

= v(rf (Q ® H T H) V (t) 

= v(r) T ((I P - lr T ) T Q(I p - lr T ) ® H T H)t](t) 

which implies, since Q is positive definite, that either ((I p — lr T ) <E) I n )v( T ) = 
or (I p ® H)t)(t) = 0. Suppose now that 

((I p - lr T ) ® Z„)r,(r) ^ . (12) 

Continuity of ?y(-) implies that there exists 5 > such that ((I p — lr T )®I n )ri(t) ^ 
for f G [r, r + §]. Therefore we must have (I p <£> H)i](t) = for t G [r, r + 6]. 
However, observability of pair (H, S) stipulates that rj(t) = for t G [r, r + <5] 
which contradicts (JH]). We then deduce ((I p - lr T ) ® I n )r](t) = for all t > 0. 
Therefore W C G ft : uj = (lr T <g> 7„)cj} = X. 

Let us now be given any solution x(-) of ([9]). Since ( that we used to 
construct Jl was arbitrary, without loss of generality, we can take x(0) = (. 
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That x(0) <E f2 implies that x(£) approaches X as t — ■> oo. Therefore we are 
allowed to write 

= lim (x(i) - (lr T ® J„)x(f)) 

t— »oo 

= lim (x(i) - (lr T ® e st )x(0)) 
where we used (fTTj) . ■ 
The following result (cf. [7]) comes as a byproduct of Theorem [T] 
Corollary 1 Consider coupled harmonic oscillators (in M. 2 ) described by 

•Ei JJi 

V 

Vi = + ^2 Hj \Vj ~ Vi) , i = 1, 2, . . . , p . 

3=1 

Oscillators synchronize for all connected T. 

6 Solution to Problem I 

In this section we use Theorem [1] in order to reach our objectives stated in 
Section [3] We first give the following fact. 

Fact 1 Let F G ]R nx ™ be a neutrally- stable matrix with all its eigenvalues re- 
siding on the imaginary axis. Then 

P := lim t- 1 [ e FTT e FT dr (13) 

is well-defined and symmetric positive definite. It also satisfies PF + F T P = 0. 

Proof. Matrix F is similar to a skew-symmetric matrix. Therefore e Ft is 
(almost) periodic [21] . Periodicity directly yields that limit in (fT3|) exists, that 
is, P is well-defined. Similarity to a skew-symmetric matrix also brings that 
inftgR |e Ft | > and sup teH |e Ft | < oo. Same goes for F T . Therefore there exist 
scalars a, b > such that al n < e F t e Ft < bl n for all t £ M. We can then write 

al n < t^ 1 / e FTr e FT dT < bl n 
Jo 
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for all t > 0. Therefore P is positive definite. Symmetricity of P comes by 
construction. Finally, observe that 



\PF + F J P\ 



= lim r 1 

t — »oo 

= lim r 1 

t — »oo 



f [e FTr e F ^F + F T e FTr e Fr )dr 



F 1 tFt 

e e 



< lim t 

t — *oo 

< lim r x (6 + l) 

t— >oo 



e FT °e F0 



whence the result follows. ■ 

Algorithm 1 Given A £ R" xn that is neutrally stable and C £ W nxn , we ob- 
tain L £ jj™ xm ag j n ows ^ m < n be the number of eigenvalues of A that 
reside on the imaginary axis. Let n-i := n — n\. If n\ = 0, then let L := 0; else 
construct L according to the following steps. 



F 
G 



Step 1: Choose U £ U nxni and W £ M nx ™ 2 satisfying 

[U W]~ 1 A[U W] = 

where all the eigenvalues of F £ R"i xn i have zero real parts. 

Step 2: Obtain P £ K™i x "i from F by (JTHJ. 

Step 3: Finally let L := UP- X {CU) T . 
Below is our solution to Problem I. 

Theorem 2 Consider systems Let Ui — Lzi where L £ W lXm is con- 

structed according to Algorithm^ and Zi is as in ([4]). Then for all network 
topologies described by connected V, solutions Xi(-) for i — 1, 2, . . . , p synchro- 
nize to 



x(t) := (r 



T ^ e A t) 



xi (0) 



where r £ R p is such that r T T = and r T l = 1. 

Proof. Let the variables that are not introduced here be defined as in Al- 
gorithm [TJ Let H := CUP- 1 / 2 and S := p^/^pp- 1 / 2 . Then (if, S) is ob- 
servable for (C, A) is detectable. Also, note that S is skew-symmetric due to 
PF + F T P = 0. 
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We let W eM" lX " and Wt e R" 2X ™ be such that 

= [U W]' 1 . 



Note then that WU = I ni) = I n2 , WW = 0, and W^U = 0. Since 

it, = Lzi, we can combine © and (TJ| to obtain 



Axi + LC 7ij (aij - 
3=1 



Let now & S K ni and r\i £ 



" 6 " 




■ pi/a o 












J„ 2 _ 
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Combining (TT4"]) and (JTSJ) we can write 

p p 

3=1 3=1 

Let L be connected and r € l p be such that r T T = 0. Then define : M>c 



11 as u>j(t) 



for i = 1, 2, . . . , p. Let w := [wj^ u). 



j , ,T 



,TlT 



(14) 
(15) 

(16) 

(17) 

o — ► 
and 



v ;= [vT V2 ■ • ■ VpV ■ Starting from (fT5|) and (jTT)) we can write 

w(t) = {T(g,e' St H T He St )w(t) + (F ® e - st H T CWe Gt )w{Q) . 

Thence 



w(t) = 0)w(0) 
where 



/ t)(L (g) e - ST iJ T CWe GT )dT 
Jo 

t) := exp ^ (L <g> e - Sa H T He Sa ) daj 



v(0) (18) 



is the state transition matrix [3], From Theorem[T]we can deduce that $(t, r) is 
uniformly bounded for all t and r. Also, for any fixed r we have lim^oo $(f, t) = 
lr T 7 ni . Moreover, e St is uniformly bounded for all t, and e Gt decays expo- 
nentially as t — > co for G is Hurwitz. Therefore we can write 

lim / <5>{t,r){T ®e~ ST H T CWe GT )dT 

Jo 

= [ ( lim $(t, r)) (F <g> e- Sr H T CWe GT )dT 

Jo Vt^oo y 

/>OC 

= / (lr T (g) 7 ni )(r ®e- ST H T CWe GT )dr 
Jo 
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Then, by (fT5|) . we can write 

lira w(t) = {lr T ® 7 m )w(0) . 

t — >oc 

Therefore solutions ^i(-) synchronize to (r T ®e )w(0). Moreover, lirut_>oo v(i) = 
for G is Hurwitz. Hence we can say that solutions synchronize to 

(r T ® e Gt )v(0). As a result, solutions Xj(-) synchronize to 

zi(0) " 
x p (0) _ 



= (r T ® e At ) 

Hence the result. ■ 



UP 



-1/2 



W 





Gt 



pl/2 C /t 



H (0) 
sp(0) 



7 Solution to Problem II 



This section, in which we provide a solution to Problem II, follows closely the 
previous one. We begin with the following algorithm. 

Algorithm 2 Given A G M. nxn that is neutrally stable and B G M. nxm , we ob- 
tain K G U™ IX ™ as follows. Let n\<n be the number of eigenvalues of A that 
reside on the imaginary axis. Let n2 '■= n — n\. Ifn\ = 0, then let K := 0; else 
construct K according to the following steps. 



Step 1: Choose U G R nxni and W G R" x ™ 2 satisfying 
[U W^Ap W] = 



F 
G 



where all the eigenvalues of F G R" lXni have zero real parts. Let U' G R" lXn 
and £R" 2X " be such that 



= [u w}- 1 



Step 2: Obtain P G K"i x "i /rom F 6y (jTSJ). 

Step 3: find/)/ Ze* X := (WB) T PW. 
Below is our solution to Problem II. 
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Theorem 3 Consider systems ([5]). Let m = Kzi where K <E j^™><™ j s con _ 
structed according to Algorithm^ and Z{ is as in ([6]). Then for all network 
topologies described by connected T, solutions Xi(-) for i — 1, 2, . . . , p synchro- 
nize to 

"a:i(0) 
x(t) := (r T ® e At ) \ 

_x p (0) 

where r S K p is such that r T T = and r T l = 1. 

Proof. Let the variables that are not introduced here be defined as in Al- 
gorithm H Let H := {P 1 / 2 WB) T and S := p^/^FP- 1 / 2 . Then {S, H T ) is 
controllable for (A, B) is stabilizable. Also, note that 5* is skew-symmetric due 
to PF + F T P = 0. 

Since ui = Kzi, we can combine ([3]) and (O to obtain 



Axi + BK ^2 7^ (ajj - a;,) 
3=1 



Let now & e K" 1 and r/ t e M™ 2 be 







■ pl/2 - 












7„ 2 _ 









Combining (jTUJ) and (|2"0)) we can write 



r]i = Gru + W^BH^ijitj-ti) 



(19) 

(20) 

(21) 
(22) 



Looking at (f2"Tj) . by Theorem [TJ we assert that solutions synchronize to 



(r J ' ® e 5t ) 



&(o) 



Now observe that — ^0 exponentially as i — > 00 for all (i, j) pairs. 

Also recall that G is Hurwitz. From (f2"2")l we can therefore deduce by input-to- 
state stability (ISS) arguments [9] that — ► as t — ► 00 for i = 1, 2, . . . , p. 
The remainder of the proof is same as that of proof of Theorem O I 
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8 Conclusion 



Let us now briefly discuss the generality of the assumptions in the paper. For 
linear time-invariant case with identical individual system dynamics, it should 
be evident that detectability (stabilizability) assumption is indispensable for 
synchronization. Regarding the neutral stability condition, it would be of great 
interest to study the synchronization of unstable systems. However, when neu- 
tral stability assumption on individual systems is relinquished, mere connected- 
ness of the network should generally not be sufficient for individual systems to 
synchronize. The reason is that, due to unstable dynamics, the trajectories will 
tend to drift apart from each other when there is no (or very little) coupling. 
The coupling strength therefore should be above some threshold to overcome 
that tendency, which requires a stronger (than connectedness) condition on the 
network topology. 



References 

[1] D. Angeli and P.-A. Bliman. Stability of leaderless discrete-time multi- 
agent systems. Mathematics of Control, Signals & Systems, 18:293-322, 
2006. 

[2] P.J. Antsaklis and A.N. Michel. Linear Systems. McGraw-Hill, 1997. 

[3] M. Arcak. Passivity as a design tool for group coordination. IEEE Trans- 
actions on Automatic Control, 52:1380-1390, 2007. 

[4] I. Belykh, V. Belykh, and M. Hasler. Generalized connection graph method 
for synchronization in asymmetrical networks. Physica D, 224:42-51, 2006. 

[5] C. Cai and G. Chen. Synchronization of complex dynamical networks by 
the incremental ISS approach. Phyisica A, 371:754-766, 2006. 

[6] Y. Chen, G. Rangarajan, and M. Ding. General stability analysis of 
synchronized dynamics in coupled systems. Physical Review Letter E, 
67:026209, 2003. 

[7] Q. Hui, W.M. Haddad, and S.P. Bhat. Finite-time semistability theory 
with applications to consensus protocols in dynamical networks. In Proc. 
of the 2007 American Control Conference, pages 2411-2416, 2007. 

[8] A. Jadbabaie, J. Lin, and A.S. Morse. Coordination of groups of mobile 
autonomous agents using nearest neighbor rules. IEEE Transactions on 
Automatic Control, 48:988-1001, 2003. 

[9] H.K. Khalil. Nonlinear Systems. Prentice Hall, 1996. 

[10] L. Morcau. Stability of multi-agent systems with time-dependent com- 
munication links. IEEE Transactions on Automatic Control, 50:169-182, 
2005. 



14 



[11] R. Olfati-Sabcr, J. A. Fax, and R.M. Murray. Consensus and cooperation 
in networked multi-agent systems. Proceedings of the IEEE, 95:215-233, 
2007. 

[12] R. Olfati-Saber and R.M. Murray. Consensus problems in networks of 
agents with switching topology and time-delays. IEEE Transactions on 
Automatic Control, 49(9):1520-1533, 2004. 

[13] L.M. Pecora and T.L. Carroll. Master stability functions for synchronized 
coupled systems. Physical Review Letters, 80:2109-2112, 1998. 

[14] A. Pogromsky and H. Nijmcijcr. Cooperative oscillatory behavior of mu- 
tually coupled dynamical systems. IEEE Transactions on Circuits and 
Systems-I, 48:152-162, 2001. 

[15] W. Ren and R.W. Beard. Consensus seeking in multiagent systems un- 
der dynamically changing interaction topologies. IEEE Transactions on 
Automatic Control, 50:655-661, 2005. 

[16] J.J.E. Slotinc, W. Wang, and K. El-Rifai. Contraction analysis of syn- 
chronization in networks of nonlinearly coupled oscillators. In Proc. of the 
16th International Symposium on Mathematical Theory of Networks and 
Systems, Belgium, 2004. 

[17] G.-B. Stan and R. Sepulchre. Analysis of interconnected oscillators by 
dissipativity theory. IEEE Transactions on Automatic Control, 52:256- 
270, 2007. 

[18] S.H. Strogatz. Exploring complex networks. Nature, 410:268-276, 2001. 

[19] J.N. Tsitsiklis, D.P. Bertsekas, and M. Athans. Distributed asynchronous 
deterministic and stochastic gradient optimization algorithms. IEEE 
Transactions on Automatic Control, 31:803-812, 1986. 

[20] S.E. Tuna. Synchronizing discrete-time neutrally stable linear systems via 
partial-state coupling. arXiv : 0801 . 2945v2 [Math. DS] , 2008. 

[21] F.S. Van Vlcck. A note on the relation between periodic and orthogo- 
nal fundamental solutions of linear systems. The American Mathematical 
Monthly, 71:406-408, 1964. 

[22] X.F. Wang. Complex networks: topology, dynamics and synchronization. 
International Journal of Bifurcation and Chaos, 12:885-916, 2002. 

[23] C.W. Wu. Synchronization in arrays of coupled nonlinear systems: passiv- 
ity, circle criterion, and observer design. IEEE Transactions on Circuits 
and Systems I, 48:1257-1261, 2001. 

[24] C.W. Wu and L.O. Chua. Synchronization in an array of linearly cou- 
pled dynamical systems. IEEE Transactions on Circuits and Systems-I, 
42:430447, 1995. 



15 



